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In this paper, we propose a disordered heterostructure in which the distribution of refractive index 
of one of its constituents follows a Levy-type distribution characterized by the exponent a. For the 
normal and oblique incidences, the effect of a variation on the localization length is investigated in 
different frequency ranges. As a result, the controllability of Anderson localization can be achieved 
by changing the exponent a in the disordered structure having heavy tailed distribution. 


I. INTRODUCTION 

In the past few decades, localization of waves in dis¬ 
ordered structures known as Anderson localization has 
attracted growing attention, in various fields of physics. 
Starting from early studies of the localization of elec¬ 
tronic wave functions in disordered crystals [I|, it is be¬ 
lieved that the localized states appear in a wide variety 
of classical and quantum materials. The possible occur¬ 
rence of Anderson localization for electrons in disordered 
solids 0, ultrasound and acoustic waves 0, the trans¬ 
port of light [1,0, microwave Q, matter waves 0, 0 
and cold atoms [9| are just a few examples among many 
others. 

All these disordered systems are typically composed 
of regions with very different kinds of the spatial distri¬ 
bution of disorder. The question is what one would ex¬ 
pect for effects of the distribution of disorder on control¬ 
ling the localization properties. This interesting question 
has been basis of numerous theoretical and experimental 
studies [IMliiMl,!!!. 

More recently, some engineered tunable random medi¬ 
ums have been assembled in the lab that follows the Levy 
type distribution [S [Bill- Levy processes have a 
power-law distribution p(x] ^ where a is the so 

called Levy exponent |l9l ]. Power law distributions have 
been appeared in different physical phenomena such as 
spectral fluctuation in random lasers [^ . [B | , superdif- 
fusive transport of light across glass microspheres whose 
diameters have levy distribution (^ . quantum coherent 
transport of electrons in one- dimensional (ID) disor¬ 
dered wires with Levy-type distribution [B| . 

Recently, Fernandez-Marin et al. have investigated 
numerically how transmission of electromagnetic waves 
varies with the system size in a ID random system in 
which the layer thicknesses follow a Levy-type distri¬ 
bution with the exponent a [l^ . They demonstrated 
that (—InT) oc for 0 < a < 1 , whereas (—InT) 
is proportional to the system size L for 1 < a < 2 
[l6|. Furthermore, Fernandez-Marn et al. have exper¬ 
imentally studied the microwave electromagnetic wave 
transmission through a waveguide composed of dielec- 
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trie slabs where spacing between them follows a distribu¬ 
tion with a power-law tail (a-stable Levy distribution) 
[B- They observed an anomalous localization for the 
case 0 < a < 1 in which transmission (T) decays with 
length of waveguide as L~°‘ while in the case of standard 
localization (T) oc exp(—L/^), ^ being the localization 
length (B| . 

More recently, Zakeri et al. have investigated Ander¬ 
son localization of the classical lattice waves in a chain 
with mass impurities distributed randomly through a 
power-law relation where s is the distance be¬ 

tween two successive Impurities [l^. For this ID har¬ 
monic disordered lattice of n-sites with random masses 
TO„, they have indicated that in the small frequencies 
for 1 < a < 2, the localization length behaves as 
C(w) -- a;-“. 

In this paper, we propose a ID disordered multilay¬ 
ered structure in which the random refractive index of 
one of its constituents follows a probability density func¬ 
tion with a power-law tail with exponent a (Levy-type 
distribution). The propagation of an obliquely incident 
electromagnetic wave through the structure is studied us¬ 
ing the transfer matrix method. The effects of variation 
of the exponent a on the localization length and Ander¬ 
son localization are investigated in different frequency 
ranges. It is shown that the localization length decreases 
with decreasing a for both 1 < a < 2 and 0 < a < 1 
at intermediate and large frequencies. Furthermore, at 
small frequencies the localization length depends on val¬ 
ues of a. In this structure, the localization length can be 
also affected by the incident angle. 

The paper is organized as follows: in Sec. II we intro¬ 
duce the notation. The numerical results are discussed 
in Sec. III. Finally, we draw our conclusions in Sec. IV. 

II. DEFINITIONS AND SETTINGS 

The disordered multilayered structure that we shall to 
study composed of an alternating sequence of layers of A 
and B having thickness of dA an dg. Fig. (HD displays a 
scheme of the proposed structure. It is assumed that the 
layers B have the same refractive index n b and the same 
width of ds, while layers A have the same thickness of 
dA and random refractive indices nAj{j = 1,2,...,V). 
The number of layers in the structure is taken to be 2N 
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FIG. 1: (Color online)Sketch of a heterostructure of length L 
with N layers with the fixed thickness dA{B)- The refractive 
index ns is fixed and the distribution of the riAS follows a 
Levy-type distribution. 


and the 2 -axis is directed across the layers. 

Here we consider UAj oc where 5j are indepen¬ 

dent identically distributed random variables with stan¬ 
dard symmetric a-stable distribution. The procedure of 
computer simulating realizations of the random variables 
5j is the following [13 ■ 


^ sin(aH) ^cos (H — oF)^ 

^ ^ [cos(H)]^/“ V ^ J ^ ^ 

where random variable W has exponential distribu¬ 
tion with mean 1 and V is uniformly distributed on 
(—7r/2,7r/2). The algorithm introduced in Eq. ([T]) al¬ 
lows us to generate a sequence of random numbers with 
a-stable distribution for the whole range of parameter 
0 < a ^ 2. The main feature of a a-stable Levy density 
distribution P(S) is the power-law, decay of its tail which 
behaves as: P(S) ~ 

We choose random numbers Sj which their absolute 
values are in the range 5 <| Sj \< 100. The refractive 
index of layer A is taken to be haj = y/| 6j \ /5. 

The transfer matrix formalism is used to compute 
the localization length of the structure. We consider a 
monochromatic electromagnetic wave obliquely incident 
from left into the random structure. In Fig. m, k and 9 
denote the wave vector and incident angle, respectively. 
The wave vector k is taken to be in xz-plane. The elec¬ 
tric and magnetic fields at incident and exit ends of the 
structure can be related by the product of transfer matrix 
of different layers included in the heterostructure as: 

= M1M2M3 ... A'hN = M (^2^+1) ,(2) 

where M is the total transfer matrix of the system and 
Mj{j = 1,2, ...,2A^) is the transfer matrix of the 


dielectric layer which is defined as: 


Mj = 


di) 


cos (^kzjdj^ sioi^k^j^^j 

—i sin (kzjdj) cos (fc 






( 3 ) 


here TyJ = {cq/ p,Q)ej/ cos'^ 9j for TM case and rjj = 
{eo/fio)ej cos^ 6j for TE case, k^j = kj cos 9j denotes the 
2 component of the wave vector in dielectric layers, and dj 
is the thickness of different dielectric layers. For a plane 
wave strikes from left into the disordered structure, the 
transmission coefficient t{uj) is expressed in terms of the 
matrix element of M as follows: 


t{uj) = 


‘^Vo 


Miirjo -I- Mi2»7oW-i-i + ^21 + M22r]N+i 


( 4 ) 


The corresponding transmittance of the structure at fre¬ 
quency UJ is: 


T(c^) = M±1 I |2 . (5) 

Vo 

To study the localization behavior in the ID random sys¬ 
tem, it is required to evaluate the localization length. 
Since the transmittance in the localized regime exponen¬ 
tially decays with the system length L, the localization 
length ^ can be numerically calculated as 


^(ta) = — lim 


2L 


L -^00 In (T(ui, L)) 


( 6 ) 


For a sufficiently long random-layered system, ^ obtained 
from the above equation is a nonrandom number due 
to self-averaging. However, for a system with a hnite 
size, the localization length can be obtained by ensemble 
averaging of the transmittance T over many realizations. 
This means that we introduce the localization length of 
a finite random configuration as 


C(w) 


2L 

(ln(r(u;,L)))- 


( 7 ) 


where (...) stands for the ensemble averaging. The val¬ 
ues of those parameters used in the following calculations 
are dA = 20mm, ds = 40mm, ns = 1 and N = 500. In 
the next section, we will discuss our numerical findings. 


III. RESULTS AND DISCUSSION 

We first study the case at which the plane wave is 
normally incident into the random structure. It is as¬ 
sumed the frequency of the incident wave is in the range 
3 X 10®(rad/s) < w < 3 x 10^°( rad/s). In order to un¬ 
derstand how localization length can be affected by the 
variation of a values, in Fig. we plot the localization 
length in units of the system length (normalized local¬ 
ization length) as a function of frequency w for different 
values of a. To obtain localization length, 2 x 10^ differ¬ 
ent random realizations with the same a values and the 
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FIG. 2: (Color online) Localization lengths in units of the 
system size for normal incidence in the frequency range 3 x 
10®(rad/s) < a; < 3 x 10^°(rad/s) for different a values. 


same number of layers are considered. As shown in Fig. 
m, the localization length decreases with decreasing the 
exponent a at some frequency ranges. For frequencies at 
which normalized localization length of the wave is lower 
than one, the system is in the localized regime. When 
a decreases from 1.8 to 0.6, the minimum frequency at 
which localization occurs shifts toward lower frequencies. 
Our numerical results confirm the self-averaging of the 
localization length. That is the localization length ob¬ 
tained from Eq. (7) does not significantly differ from 
the localization length of a single realization with large 
number of layers. Furthermore, increasing the number of 
random realizations from 2 x 10^ does not lead to any 
considerable change in localization length values. 

To determine whether or not this dependence of local¬ 
ization length on the a values is observed in the case of 
oblique incidence, we display the normalized localization 
length versus frequency for Q = 15° and 6 = 30° in Fig. 

The corresponding frequency range is the same as 
that in Fig. ([2]). It is clearly seen in Fig. m that for 
oblique incidence the localization length decreases when 
the exponent a decreases from 1.8 to 0.6. To better un¬ 
derstand the effect of incident angle on the localization 
length, in Fig. we show the normalized localization 
length versus a for different incident angles 9 = 0, 15° 
and 30° with the same a = 1.8. 

One can see that the localization length increases with 
incident angle at some frequency regions. Furthermore, 
the dip in localization length moves to higher frequencies 
with increasing incident angle. Therefore, our calculated 
results indicate that the localization length depends on 
the incident angle as well as the exponent a. 

It is well known that for TM polarized waves propagat¬ 
ing in a ID random structure the localization length takes 
a large maximum value at some critical angles which 
are called generalized Brewster angles [2^. It has been 
demonstrated that generalized Brewster angle increases 
from 0° to 90° with increasing the disorder-averaged re¬ 
fractive index [2^. The localization length for a weak 
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FIG. 3: (Color online) Localization lengths in units of the 
system size for oblique incidence 9 = 15° (top) and 9 = 30° 
(bottom) in the frequency range 3 x 10®(rad/s) < ca < 3 x 
10^° (rad/s) for different a values. 


disorder diverges in the vicinity oh 0 = 45° when the 
average of refractive index is equal to 1 [H, ■ This 

phenomenon is known as Brewster anomalies and the cor¬ 
responding angle is called the Brewster angle [2^, . In 

Fig. ®, we display the localization length versus inci¬ 
dent angle for a = 0.6 and a; = 2 x 10®(rad/s). One 
can see that the generalized Brewster angle is about 80° 
at which the localization length is significantly enhanced 
over its value at the normal incidence {9 = 0°). More¬ 
over, for incident angle in the range 68° < 9 < 88°, the 
system is in the extended regime, while at other incident 
angle the system is in the localized regime. The effects of 
a variation on the generalized Brewster angle and Brew¬ 
ster anomalies are under investigation and their results 
will be reported in the near future. 

To understand the physical reason for the behavior of 
localization length with a, in Fig. ®, we plot the mean 
value of refractive index versus a. As shown in Fig. ®, 
the average value of refractive index increases with in¬ 
creasing a. As a result, increases of a leads to the de¬ 
creasing the refractive index contrast between the layers 
of disordered system. This effect causes the scattering 
strength to decrease. Hence, we expect that the localiza¬ 
tion length decreases with increasing a. In addition, the 
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FIG. 4: (Color online)Localization lengths in units of the sys¬ 
tem size for different incident angles 0 — 0, 15° and 30° in the 
frequency range 3 x 10® (rad/s) < w < 3 x 10^° (rad/s) with 
the same a values. 



FIG. 5: (Golor online)Localization lengths in units of the 
system size versus incident angle 9 for a = 0.6 and oj = 
2 X 10®(rad/s). 


Fig. dH) represents the variance of refractive index ver¬ 
sus a. One can see that the variance of refractive index 
increases with decreasing a. Therefore, the randomness 
strength increases with decreasing a. This effect also re¬ 
sults in the enhancement of localization with decreasing 
a. 

Now, we shall to study how heavy-tail distribution of 
random refractive index ua affects the normalized local¬ 
ization length in the frequency range 2.5 x 10^^(rad/s) < 
w < 3 X 10^^(rad/s). Fig. d?)) shows the corresponding 
results for different values of a = 1.8, 1.4, 1.0 and 0.6. 

It is clearly seen that the normalized localization 
length shows an oscillatory behavior in this frequency 
range for different a values. Moreover, for all frequencies 
in this frequency range and for all a values, we have a 
localized mode whose localization length decreases with 
decreasing a. As a result, decreasing a value improves 
the localization. This effect is attributed to the enhance¬ 
ment of mean value and variance of refractive index with 


FIG. 6: (Golor online) Mean and variance (inset) values of 
refractive index ua as a function of the exponent a. 



FIG. 7: (Color online) Localization lengths in units of the 
system size for normal incidence in the frequency range 2.5 x 
10^^ (rad/s) < oj < 3 x 3^^ (rad/s) for different a values. 


decreasing a. 

Next, we consider the effect of a variation on the nor¬ 
malized localization length for longer wavelengths. The 
normalized localization lengths versus uj in frequency 
range 3 x 10®(rad/s) < a; < 2.5 x 10®(rad/s) are dis¬ 
played in Fig. dS]) for different values of a = 1.8, 1.4, 1.0 
and 0.6. 

As shown in this figure, the localization length indi¬ 
cates an oscillatory behavior where their peak values in¬ 
creases with lowering the frequency for all a values. For 
all wavelengths, the system is in the extended regime. 
With decreasing the a values, the peaks of normalized 
localization length shift to lower frequencies, while the 
peak values decreases. 

Consequently, our results demonstrate that if the dis¬ 
tribution of the random layer in a ID disordered system 
follows the a-stable Levy distribution, the localization 
length can be manipulated with the exponent a. Due 
to the random nature of Anderson localization, the con¬ 
trol of this phenomenon in a regular manner is of great 
importance and has potential applications. 
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FIG. 8: (Color online) Localization lengths in units of the 
system size for normal incidence in the frequency range 3 x 
10®(rad/s) < u < 2.5 x 10®(rad/s) for different a values. 

IV. CONCLUSION 

We have studied the localization of an electromagnetic 
wave normally and obliquely incident into a ID disor¬ 
dered structure where the refractive index of one its con¬ 
stituents is fixed while that of the other constituents 
is a random number drawn from a Levy-type distribu¬ 
tion with exponent a. It has been demonstrated that 


for normal incidence the localization behavior in this 
structure can be manipulated in a regular manner by 
changing the value of a. When a decreases, the local¬ 
ization length of the waves decreases in the frequency 
range 3 x 10®(rad/s) < a; < 1.9 x 10^°(rad/s) decreases. 
This effect is due to the increase of mean and variance of 
refractive index with decrease of a. Moreover, the min¬ 
imum frequency at which localization appears shifts to 
lower frequencies with decreasing a. Localization length 
shows the same trend with variation of a for oblique in¬ 
cidence. It has been also investigated how the localiza¬ 
tion length can be affected by a variation in the lower 
and higher frequency ranges. In the frequency range 
2.5 X 10^^(rad/s) < a; < 3 x 10^^(rad/s), the system is in 
the localized regime and the localization length indicate 
an oscillatory behavior with approximately fixed ampli¬ 
tudes. Decrease of a in this frequency range also gives 
rise to the enhancement of localization. In the frequency 
range 3 x 10®(rad/s) < w < 2.5 x 10®(rad/s). it is found 
the system is in the extended regime and the localization 
length exhibits an oscillatory behavior with increasing 
amplitude whose value decreases with decreasing a. Fur¬ 
thermore, reduction of a causes the peaks of localization 
length to shift toward lower frequencies. Consequently, 
in disordered media, employing the a-stable Levy distri¬ 
bution provides a way to easily control the localization 
phenomenon. 
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